Universality of massive scalar field late-time tails in black-hole spacetimes 
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The late-time tails of a massive scalar field in the spacetime of black holes are studied numerically. 
Previous analytical results for a Schwarzschild black hole are confirmed: The late-time behavior of 
the field as recorded by a static observer is given by ip{t) ~ t~ 5 ^ 6 sin[u>(t) x t], where u>(t) depends 
weakly on time. This result is carried over to the case of a Kerr black hole. In particular, it is found 
that the power-law index of —5/6 depends on neither the multipole mode £ nor on the spin rate 
of the black hole a/M. In all black hole spacetimes, massive scalar fields have the same late-time 
behavior irrespective of their initial data (i.e., angular distribution). Their late-time behavior is 
universal. 
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I. INTRODUCTION 

The late-time behavior of massless fields in black hole spacetime has been studied in detail for both linearized 0, Q 
and fully nonlinear (spherical) evolutions 0,0. In contrast, the late-time behavior of massive fields has been studied 
in much less detail. In this paper we study the late-time behavior of a scalar field numerically, in the spacetimes of 
Schwarzschild and Kerr black holes. 

The first to consider the problem of a massive scalar field in the spacetime of a Schwarzschild black hole were 
Novikov and Starobinski, who studied the problem in the frequency domain, and found that there are poles in the 
complex plane closer to the real axis than in the massless case. They thus inferred that the decay rate would be slower 
in the massive case than in the massless case [EJ. That problem was later studied numerically by Burko both for 
a linearized massive scalar field, and for a self-gravitating, spherical massive field. In both cases it was found in [(| 
that the late-time behavior was given by an oscillating field, whose envelope decayed according to an inverse power 
law, and whose frequency was determined by the Compton wavelength of the field, i.e., by its mass term. In [™it 
was also reported that the decay rate of the envelope was given by t~ a , where a ~ 0.8. No attempt was done in |(J, 
however, to determine the value of a accurately, or to determine its dependence or lack thereof of the parameters of 
the problem. The most detailed analytical study of the problem was done by Koyama and Tomimatsu in a series of 
three papers 0,H@. 

The most striking feature of the problem of tails of a massive scalar field is that the tails exist already in flat 
spacetime |ll|. because spacetime acts like a dispersive medium for a Klein-Gordon field. Based on the exact Green 
function, the behavior of the field in Minkowski spacetime was found in jfj to be given by 

V'flat - i-'" 3/2 sin(^) , (1.1) 

£ being the multipole moment of the field, and being its Compton wavelength. A derivation of Eq. I|l.l|) appears 
in the Appendix. More recently, an interesting self similar behavior of a Klein-Gordon field in Minkowski spacetime 
was found in [lj} ■ The self similar behavior reduces to Eq. i|l.lfl for a field evaluated at a fixed location at very late 
times. 

In the spacetime of a Schwarzschild black hole there are important differences. The decay rate was found in |8( to 
be given by 

V'sch ~ i~ 5/6 sin[w(t) x t] , (1.2) 

where the power index is independent of I. The time-evolving angular frequency u>{t) approaches /j, asymptotically as 
t — ► oo, but is different from fi at finite values of time. This result was generalized in |9| to any spherically symmetric, 
static black hole spacetime. The details of this result have not been corroborated numerically — although qualitatively 
they are supported by the results of j(| — a task we undertake in this paper. In addition, we also show that this result 
remains correct also for a Kerr black hole, for which the power index —5/6 is independent of both £ and the spin rate 
of the black hole a/M. We also study the behavior of oj(t), and find it to be in agreement with the results of H, ||. 
Similar conclusions were obtained recently also for a spinning dilaton black hole ^2 (see also fifty). 
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II. MASSIVE TAILS IN SCHWARZSCHILD SPACETIME 



Massive fields in Schwarzschild spacetime were studied numerically (for both a linearized and for a spherical self- 
gravitating Klein-Gordon field) in [g and analytically in As discussed in for fj,M < 1, the tail regime satisfies 
t ^> fi~ 3 M~ 2 , M being the mass of the black hole, and t being the regular Schwarzschild time coordinate. Therefore, 
for a light scalar field, the heavier the field is, the earlier the tails are seen. (For heavy fields, \xM ^> 1, the tail regime 
is for t ^> M, independently of the mass of the field.) To facilitate numerical simulations, one is tempted then to 
increase the mass of the field. However, increasing the mass of the scalar field requires also to increase the resolution, 
as the wavelength of the field is inversely proportional to its mass, and as one needs at least a few data points per 
wavelength to resolve an oscillating field. One faces then a competing effect: to decrease the computational time to 
an optimum one needs to balance the scalar field mass and the resolution of the numerical code. (We do not discuss 
here the case of /iilf 3> 1 as the required numerical resolution is beyond current practical computational limits.) One 
may gain insight into the above discussion, noticing, following that for // <C 1/M, the corresponding Compton 
wavelength of the scalar field A(= M x ) ^ M, that is the field's wavelength is much longer than the typical radius 
of curvature of spacetime, or the scale of inhomogeneity of curvature. We therefore expect the field at early times to 
evolve similarly to its evolution in flat spacetime. At later times the curvature effects become apparent, and we expect 
deviations from the flat-spacetime behavior. At asymptotically late times we expect to find the behavior predicted in 
The greater the mass term /i, the shorter the associated Compton wavelength A, and the sooner the asymptotic 
domain is expected. We therefore expect an oscillating field, whose envelope is a broken power law: at early times it 
is described by t~ e ~ 3 ^ 2 and at late times by i~ 5 / 6 . 

The field ip satisfies the Klein-Gordon equation 



(□-^ = 0, 



(2.1) 



□ being the DAlembertian operator in curved spacetime. We define ip = rip in the usual way, and decompose the 
field if> into Legendre modes ip = J2e iptPe{cos6). The radial equation that each mode ipt satisfies is then given by 



1 



0. 



(2.2) 



where the effective potential 



Vt(r) 



i(£ + 1) 



2M 

T3~ 



(2.3) 



To test our expectations in the context of a Schwarzschild black hole, 
solved the radial equation for We specified a pulse of the form 



(u - ui)(u - u 2 ) 



1 8 



(u 2 — U\) 2 



we used a double- null code in 1+1D, and 



(2.4) 



along the ingoing segment of the characteristic hypersurface v = for u\ < u < u%, and ip{u) = otherwise. We 
also took ip(v) = along u = 0. Here, u (v) is the usual retarded (advanced) time. The results reported on below 
correspond to u± = 20M and u 2 = 60M, but our results are unchanged also for other choices of the parameters, or 
for other choices of the characteristic pulse. In the following we observe the fields along r* = 0, r* being the usual 
Regge- Wheeler 'tortoise' coordinate. 

The code is stable, and exhibits second order convergence. The convergence test is presented in Fig.^ which shows 
the ratio of differences for three runs with different resolutions. 

We first chose a low value of fiM. Figure [5] shows the field along = const as a function of time. The initial field 
here is spherical (I = 0), and a value of \xM = 10 -3 was chosen. Similar results are obtained also for other low values 
of fj,M. After the prompt burst and the quasi-normal epoch, an oscillatory tail is seen. The decay rate of the envelope 
of the tail is very close to t~ 3 / 2 . 

We next consider a high value of \iM. Figure |3| shows the field for the same parameters as in Fig. El except that 
here fj,M = 1. Again, an oscillatory tail is seen, and the decay rate of the tail is very close to i~ 5 / 6 . 

The most interesting case is that of an intermediate value of for which both types of behavior coexist. In 
Figure 01 we show the field for fiM = 2.5 x 10~ 2 . This figure shows initial oscillatory decay with envelope decaying 
like i -3 / 2 , which changes gradually to an oscillatory decay with decay rate of t -5 / 6 . This is the broken power-law 
expected: at early times the field has not noticed yet the presence of the black hole, and therefore behaves as in flat 
spacetime (decay rate of i -3 / 2 ). At late times, the field takes its asymptotic decay rate of i~ 5 / 6 , and the onset of the 
asymptotic behavior can be controlled by changing the value of fiM . 
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FIG. 1: Convergence test for the 1+1D code. We find the times at which the field has zeros for three different resolutions, high 
t(h), medium t(m), and low t(\), and compute the ratio [t(h) — t(m)]/[t(m) — t(Yj] for each zero. We then plot this ratio as a 
function of time t. Here, we used N = 5, 10, and 20 grid points per M for the three resolutions. Second order convergence 
corresponds to the ratio equaling 0.25. 




FIG. 2: Field of a massive scalar field along r* = 0. The initial data are for i = 0, ui = 20M and U2 = 60M, and N = 1. The 
solid curve is the massive scalar field (with fj,M = 10 -3 ), and the dashed line is proportional to (t/M)~ 3 ^ 2 . 



The prediction of Ref. [8j is that the asymptotic decay rate of i~ 5 / 6 is independent of the value of £. We check this 
in Fig. which shows that same data as in Fig. [31 except that here we take I = 1. Again, the decay rate is i~ 5 / 6 , 
the same as with I = 0. We could not find any deviation from the i -5 / 6 behavior for any value of £. We therefore 
conclude that our simulations are in agreement with the prediction of Ref. Q . 

In flat spacetime the period of the oscillations does not change at late times, and is given by Tfl at = 2n/fi (see the 
Appendix). In Schwarzschild spacetime, the period is no longer fixed. In Ref. @, the oscillatory part of the field is 
given by 



sin 



/xt - (3/2)(2tt^M) 2/3 (^) 1/3 + smaller terms 



To interpret this prediction as a changing period, we re- write the oscillatory part of the field as sin[cij(i) x t], where 

2/3 /x 



Uj(t) = fl 



2 l t \t 



(2.5) 



Figure El shows the power spectrum of the field, for the case /j,M = 1. A sharp peak at / = (2ir) 1 is seen. This 
peak corresponds to angular frequency equaling the mass term /i. Notice the one-sided broadening of this peak. While 
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FIG. 3: Field of a massive scalar field along r* = 0. The initial data are for I = 0, Ui = 20M and U2 — 60M, and TV = 4. The 
solid curve is the massive scalar field (with jiM = 1), and the dashed line is proportional to (i/M) _5//6 . 
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FIG. 4: Field of a massive scalar field along r* = 0. The initial data are for I = 0, ui = 20M and U2 = 60M, and N = 4. 
The solid curve is the massive scalar field (with \iM = 2.5 x 10~ 2 ). The dashed line (a) is proportional to (t/M) -3 / 2 , and the 
dotted line (b) is proportional to (t/M)~ s ^ 6 . 

no frequencies higher than (27r) _1 appear to be present, lower frequencies are. This suggest that as time progresses, 
the frequency increases to its asymptotic value of (2n)~ 1 . We test this in Fig. 0by plotting, 1 — 2n/P(t) as a function 
of t -1 . Here, P(t) is the local period of the oscillations. We find that as t — > oo, indeed the period decreases to its 
asymptotic value. Next, we consider in greater detail the rate at which the period approaches its asymptotic value. 
Recall that according to Ref. 8], the slope of the curve in Fig.[7]should be 2/3. In Fig.|S]we display the local slope of 
the curve as a function of t . Our numerical result, of an asymptotic slope of 0.66 ± 0.01, is in agreement with the 
predictions of Ref. [8j . 

III. MASSIVE TAILS IN KERR SPACETIME 

Now we turn to the discussion of massive scalar field tails in Kerr background. We intuitively expect the late-time 
behavior in this case to be identical to that of the Schwarzschild case, as presented in the previous section. The reason 
why the late-time tail is expected to be independent of a/M is the following: For a massless field in Kerr spacetime, 
there is a mixing of modes and each existing mode decays with a decay rate of t~( 2e+3 \ which is the same decay rate 
as in the Schwarzschild case. All the modes which are not disallowed are excited, such that the overall decay rate is 
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FIG. 5: Same as Fig. 01 but with 1 = 1. 
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FIG. 6: Power spectrum of the field. Here, fiM = 1, and N = 100. To obtain this plot, the field's values were first scaled by 
£ 5//6 to get close to a pure sinusoid, and then the Fourier transform was taken. We assume here (and below) M = 1. 




FIG. 7: The changing period: 1 — 2tt/P(£) vs. M/t, for the same data as in Fig. [S| 



6 




0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 

M/t 

FIG. 8: The local slope of the curve in Fig. pit), as a function of M/t. 



dominated by the existing mode with the slowest decay rate @] ■ We expect this situation to remain basically the same 
also for a massive field. However, for a massive scalar field in Schwarzschild spacetime the decay rate is independent 
of the mode t. This implies that although in Kerr mode mixing will indeed generate more modes and each mode will 
have the same evolution as in Schwarzschild spacetime (as in the massless case) , because in the Schwarzschild case all 
massive field modes have the same decay rate, they will all have the same decay rate in Kerr spacetime. 

Our numerical simulations were performed using the penetrating Teukolsky code (PTC) [lj], which solves the 
Teukolsky equation for linearized perturbations over a Kerr background in the ingoing Kerr coordinates (t,r,9,tp). 
These coordinates are related to the usual Boyer-Lindquist coordinates (t,r,9,(p) through (p = cp + JaA" 1 dr and 
t = t — r + r„, where A = r 2 — 2Mr + a 2 and = J (r 2 + a 2 ) A -1 dr. The Teukolsky equation for the function ip in 
the ingoing Kerr coordinates can be obtained by implementing black hole perturbation theory (with a minor rescaling 
of the Kinnersley tetrad [lj]). It has no singularities at the event horizon, and therefore is capable of evolving data 
across it. This equation is given for a massless field (/i — 0) by 



d 2 tp d 2 tp 



dtp 



+ - A^+2(*-l)(r-M)^ 

dtdr 



dt 2 

i d ( dtp 



sin 6 d<p 



-2a 



d 2 tp 2s cot 8 dtp . , , „ 

— - i . „ ttz + (s COt 2 6 + s)lp 



drdip 
+2 [sr + ias cos 



» dtp 
(.-1)^=0. 



(3.1) 



For the case of interest, we set s = in Eq. and include a mass term /j, 2 (r 2 + a 2 cos 2 6)ip on the left hand 

side of Eq. I|3.1|l . The PTC implements the numerical integration of the resulting equation by decomposing it into 
azimuthal angular modes and evolving each such mode using a reduced 2+1 dimensional linear partial differential 
equation. The results obtained from this code are independent of the choice of boundary conditions, because the inner 
boundary is typically placed inside the horizon, whereas the outer boundary is placed far enough that it has no effect 
on the evolution. (As was shown in Ref. |l5j . close timelike boundaries with Sommerfeld-like boundary conditions 
do not allow for the evolution of late-time tails.) Typically, for the simulations performed in this study, the outer 
boundary was located at 4000M and a grid of size 40000 x 32 (radial x angle) was used. The initial data was always 
chosen to be a gaussian distribution, centered at 50M and with a width of 2M. The PTC is stable and exhibits 
second order convergence as clear from Fig. O 

We first demonstrate the independence of the late-time evolution of the Kerr parameter (a/M) of the background 
spacetime. Figure 1101 (left panels and upper panel on the right) shows tails from several different evolutions corre- 
sponding to different values of a/M, all plotted together. Each oscillatory tail shown has the expected period (about 
2n/(j.) and a decay rate close to t~ 5 ^ 6 . The value of fiM is chosen here to be 1. 

We next demonstrate the independence of the oscillatory tail of the value of I. Figure (right panels) shows the 
late-time evolution of the massive scalar field [p,M = 1) in a Kerr background spacetime with a/M — 0.6 for several 
different values of I. There appears to be no deviation from the i -5 / 6 oscillatory tail for any value of I. 
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FIG. 9: Convergence test for the 2+1 D penetrating Teukolsky code. We compute the times at which is field has zeros for 
three separate resolutions [high t(h), medium i(m), and low i(l)] and plot the ratio [i(h) — t(m)]/[t(m) — t(l)] as a function of 
time. The resolutions used for this test, were M/40, M/20 and M/10 for the high, medium and low, respectively. Second order 
convergence is clear by the value of the convergence ratio (0.25). 

Last, we turn to the changing period of the oscillations in the Kerr case. In Fig. II II we plot the power spectrum of 
the field for the case of a/M = 0.6, £ = and /J.M = 1. Much like in the Schwarzschild case in the previous Section, 
we observe a one-sided broadening of the peak at (2tt)^ 1 . To test whether this indicates a monotonic increase in 
frequency we plot in Fig. ^1 1 — 2n/P{t) as function of t -1 , where Pit) is the local period. Indeed, we see in the 
Kerr case a decreasing period, with an asymptotic rate (slope of the curve in Fig. ^] is shown in Fig. I13|) of about 
0.68 ± 0.02 which is consistent with the (Schwarzschild) value of 2/3. 

Thus, in accord with our expectations, the asymptotic late-time behavior of massive scalar fields in a Kerr spacetime 
is identical to that in Schwarzschild background. Specifically, the asymptotic late-time tail is insensitive to the 
multipole number £, and also insensitive to the spin rate of the black hole. The late-time behavior of massive scalar 
fields in black hole spacetimes is universal. There is yet to be an analytic study in the Kerr context, akin to the one 
in the Schwarzschild case Also, it should be noted that all the numerical simulations reported on in this 

work were axisymmetric. In a nonaxisymmetric evolution of massive scalar fields in a Kerr spacetime, under certain 
conditions, a very interesting instability arises that has been studied in the frequency-domain [l6j . We are hoping to 
return to an extensive study of that case in the time-domain elsewhere. Lastly, intriguing deviations from the results 
reported on here, which appear to introduce a certain low-frequency modulation at very late times, were reported on 
in numerical simulations |l7l Hs[ . Notably, the analytical analyses do not show any such phenomena. It remains an 
open question whether such phenomena are real physical ones, or a numerical artifact. 
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APPENDIX A: TAILS OF A MASSIVE SCALAR FIELD IN FLAT SPACETIME 



In this Appendix we derive the decay rate of Klein-Gordon tails in flat spacetime. 
The exact Green's function is known , and is given by 



/ S(t-r) /xJi /zVt 2 - r 2 ] 

g (r, t) = — B(t - r) , 

r vi 2 - r 2 



(Al) 



for a source at the origin (in space and time). Here, J\ is the Bessel function of the first king of order 1. This Green's 
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FIG. 10: Massive scalar field (with fiM = 1) sampled at r = 50M. The initial data are a gaussian distribution, centered at 
50M and has a width of 2M. The dotted curves show the behavior of the scalar field, and the solid lines are proportional to 
(t/M)- 5/6 . The sequence of plots demonstrate a clear independence of a/M, and also of the mode number I. 



9 




f M 



FIG. 11: Power spectrum of the field evolving in a Kerr spacetime. Here, [t,M = 1, and a/M = 0.6. To obtain this plot, the 
field's values were first scaled by t 5 ^ 6 to get close to a pure sinusoid, and then the Fourier transform was taken. 
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FIG. 12: The changing period in Kerr spacetime: 1 — 2ir/P(t) vs. M/t for the same parameters as in Fig. 1111 
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FIG. 13: The local slope of the curve in Fig. 1121 p(t), as a function of M/t. 
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FIG. 14: The construction of the "image charge" solution. The burst of the spherical wave is at the event (0,ro). The domain 
of dependence of this event is the union of regions I and III in the diagram. As this solution is not regular along r — we 
introduce an "image charge," i.e., a burst with the opposite sign at the fictitious event (0, — ro). The field due to this image 
solution in region II is of course unphysical, but in region III it is, and superposes with the solution of the original burst in 
that region. 



function can of course be integrated to find the field straightforwardly, but instead we shall use a more geometric 
approach to find the field. 

Consider a spherical shell of radius ro, such that the initial perturbation field is well localized in space and in time. 
The only requirement is that the perturbation will be regular everywhere, and in particular at the location and time 
of the burst and at the origin r = 0. The field <f> satisfies the Klein-Gordon equation 

(a-fi 2 )4> = 0. (A2) 

We next define, as usual, ip = r( t ) : an d write the radial equation for ip. 

As the field propagates in spherical waves, the field can depend only on the spacetime interval s from the event of 
the burst. That is, ip — ip(s). Inserting this into the radial equation, we find that ip satisfies 

s V + sip' + ii 2 s 2 i> = , (A3) 

a prime denoting derivative with respect to the interval s. This is just the Bessel equation, with general solution 

ip(s) = ci Joi^s) + c 2 Y (fJ,s) . (A4) 

Here, Jo and Y"o are the Bessel functions of the first and second kinds, respectively, of order 0. The integration constant 
£2 must vanish for the field to be regular at the event of the burst (s = 0). Therefore, 

ip(s) = cxJoGtzs) . (A5) 

For r = const and for t 3> r, ro, we then find that 

cos( rit — 7r/4) 

uJ ■ (A6) 

This solution does not satisfy the regularity requirement, as the field <f> diverges at the spatial origin r = 0. However, 
the field <fi is required to be regular at r = 0, especially at late times (/it 3> 1). 

To obtain the physical solution we superpose solutions to the radial equation, such that the field is regular at the 
origin. This can be achieved by an "image charge." The superposition 

ip(s) = ci [Jo(M s O - Jo(ns 2 )] , (A7) 

is the requested solution, where sf = t 2 — (r — r ) 2 and s 2 , — t 2 — (r + r ) 2 . This corresponds to an "image charge" at 
r = — ro. Obviously, r = — ro is not a physical point in spacetime. It is a mathematical way to write a formal solution, 
where in the physical part of the domain of dependence of the superposed bursts the solution is the one sought (see 
Fig. I14fl . In fact, this is just a mathematical way to specify regular boundary conditions along r = 0. 
For r = const and for t 3> r, ro , we then find that 

cos(ut — 7r/4) , . . 

i>(t » r, r ) « c x // V r . (A8) 
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The phase is unimportant, such that we finally get the result 

0(t»r,r o )«cirop^ 7 |. (A9) 

The fall off rate of i~ 3 / 2 is the one reported on in Ref. [Tl|. 

To find the time dependence of higher £ modes, we find the derivative of the spherical solution with respect to the 
Cartesian coordinate z. Recall that on the equatorial plane, dPi(cos8)/ dz = — (£ + l)Pi + i(cos9)/r. We can thus 
generate ip\ from tpo by calculating dipe = o/ dz, and by the uniqueness of the solution that should give us the dipolc 
solution tpe—i. Specifically, notice that 

d dr ds d r d 
dz dz dr ds s ds 

Notice the factor 1/s in this expression. This implies that whenever a derivative with respect to z is taken, the 
exponent of t in the denominator increases by 1. (This is the case because of the sinusoidal function in the numerator. 
The leading term in i" 1 comes from its differentiation, and then the exponent of the denominator does not change. 
It is increased by 1 because of this 1/s factor.) 

Carrying i times this differentiation with respect to z, we find that 

sin(ui) 

^(*»r,r )~-£^. (A10) 

The reason why the Huygens principle fails in this case is that the phase velocity of a plane wave satisfying the 
Klein-Gordon equation is 

*"7TOT (AU) 

so that whenever /i ^ 0, different frequencies which make up a wave packet travel at different speeds. As noted in 
[Tl| , it is not surprising that the associated plane waves do not arrive at the observation point with the same relative 
phase they had at the beginning. In fact, spacetime behaves like a dispersive medium for the Klein-Gordon equation, 
as the wave number k = wyl — (h/oj) 2 is no longer a linear function of u. 
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